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Abstract —  Revisiting  an  earlier  examined  multiple-input 
multiple-output  (MIMO)  radar  waveform  design  problem  which 
optimizes  both  minimum  mean-square  error  estimation  (MMSE) 
and  mutual  information  (MI),  we  formulate  a  new  waveform 
design  problem  and  provide  some  further  results  in  this  paper, 
which  complements  the  previous  study.  More  specifically,  we 
present  an  iterative  optimization  algorithm  based  on  the  alternat¬ 
ing  projection  method,  to  determine  waveform  solutions  that  can 
simultaneously  satisfy  a  structure  constraint  and  optimize  the  de¬ 
sign  criteria.  Numerical  examples  are  provided,  which  illustrate 
the  effectiveness  of  the  proposed  approach.  In  particular,  we 
find  that  the  waveform  solutions  obtained  through  our  proposed 
algorithm  can  achieve  very  close  and  virtually  indistinguishable 
performance  from  that  predicted  in  the  previous  study. 

I.  Introduction 

As  multiple-input  multiple-output  (MIMO)  radar  has 
emerged  as  a  new  paradigm  in  the  modern  radar  technology 
(see  [1]  for  a  review),  a  lot  of  efforts  have  been  devoted 
to  designing  waveforms  for  MIMO  radar  to  unleash  its  ex¬ 
cellent  performance  advantages.  Among  the  existing  works 
is  [2],  which  studies  MIMO  radar  waveform  design  for  the 
identification  and  classification  of  extended  radar  targets.  It 
focuses  on  optimizing  two  criteria:  maximization  of  the  mutual 
information  (MI)  and  minimization  of  the  minimum  mean- 
square  error  (MMSE),  and  requires  the  a  priori  knowledge 
of  the  extended  target  impulse  response,  e.g.,  the  target’s 
second-order  statistics.  With  such  knowledge  and  under  some 
conditions,  [2]  indicates  that  under  the  same  constraint  that 
limits  the  total  power,  these  two  criteria  lead  to  the  same 
optimum  solution,  which  employs  water-filling  and  allocates 
the  transmitted  power  in  proportion  to  the  quality  of  the 
particular  mode  in  question.  In  a  sequent  paper  [3],  the 
minimax  robust  waveform  design  is  investigated,  which  leads 
to  different  optimal  waveforms  under  the  two  criteria  above. 
In  both  [2]  and  [3],  some  fundamental  quantity  has  been 
determined  to  characterize  the  essential  part  of  the  optimum 
waveform  design,  which  can  be  deemed  as  the  covariance 
matrix  of  the  transmitted  waveforms.  However,  unfortunately. 
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the  solutions  obtained  in  [2]  [3]  are  not  in  the  ultimate  form 
of  the  transmitted  waveform  matrix.  Further,  the  solutions 
identified  in  [2]  [3]  require  a  certain  condition  embodied  in,  for 
example.  Theorem  1  of  [2],  that  the  waveform  matrix  should 
have  a  specific  Kronecker  structure.  Herein  lies  the  purpose 
of  this  work:  to  find  alternative  formulations,  without  any 
restrictions,  which  provide  directly  the  transmitted  waveform 
matrix. 

By  explicitly  taking  into  account  the  conditions  required 
by  Theorem  1  of  [2],  in  this  paper  we  reformulate  the  MIMO 
radar  waveform  design  problem  that  was  studied  in  [2]-[4], 
We  further  present  an  iterative  optimization  algorithm  to  seek 
waveform  solutions  which  both  satisfy  the  structure  condition 
and  optimize  the  design  criteria.  Such  algorithm,  assuming 
the  method  of  alternating  projection  as  a  basis,  cyclically 
and  alternatively  computes  the  best  approximations  under 
different  constraints.  It  also  involves  solving  several  matrix 
nearness  problems,  for  which  we  resort  to  methods  such  as 
the  Procrustean  transformation  and  the  Kronecker  product 
approximation.  We  also  provide  some  numerical  examples, 
through  which  we  find  that  the  waveform  solutions  we  obtain 
by  applying  the  proposed  algorithm,  though  are  no  longer  in 
closed  form,  lead  to  performance  which  is  very  close  to,  and 
almost  indistinguishable  from  that  predicted  by  [2]  [3], 

The  remainder  of  this  paper  is  organized  as  follows.  In 
Section  II,  we  introduce  the  signal  model  and  summarize  the 
results  about  the  optimal  waveform  design  originated  in  [2], 
In  Section  III,  we  reformulate  the  waveform  design  problem. 
In  Section  IV,  we  present  an  iterative  optimization  algorithm 
for  the  waveform  design.  In  Section  V,  we  provide  numerical 
examples  to  validate  the  proposed  waveform  design.  This 
paper  is  concluded  in  Section  VI. 

Notation :  In  this  paper,  we  use  the  symbol  ||  •  ||^  for  the 
Frobenius  norm  of  a  matrix,  (-)+  for  the  nonnegative  part  of 
a  real-valued  matrix,  and  (g)  for  the  Kronecker  product.  ‘ft{  ■ } 
and  { ■ }  denote  the  real  part  and  the  conjugate  of  a  complex 
number,  respectively. 

II.  Signal  Model  and  Optimal  Waveform  Design 

We  consider  here  a  MIMO  radar  system  with  P  transmit 
and  Q  receive  antennas  which  may  be  distributed  in  space. 
We  model  the  reflection  from  the  signal  sent  from  the  pth 
transmit  element  and  captured  at  the  gth  receive  element  using 


a  finite  impulse  response  (FIR)  linear  system  with  order  v  in 
baseband,  whose  impulse  response  is  g<'p’q\l)1  l  =  0,  ■  •  •  ,  v. 
Then  the  received  waveform  at  the  q\h  receive  element  and 
discrete-time  k  can  be  modeled  as 

P  V 

Vq{k)  =  EE  g(p,q\l)  xp(k  -  l)  +  nq(k)  (1) 

p= 1  1=0 

where  xp(k)  is  the  waveform  transmitted  from  the  pth  trans¬ 
mit  element  and  nq(k)  is  the  additive  complex  Gaussian 
noise  measured  at  the  gth  receive  element.  Let  L  denote 
the  length  of  the  observed  signal  vector  starting  at  an  ar¬ 
bitrary  discrete  time  k,  and  we  assume  that  L  >  PM 
(typically  L  PM),  where  M  =  v  +  1.  Denote  g <'p'q'>  = 
[s(p’9)(0),  ■■■  ,  g(p'q) {v)]T ,  yq  =  [ yq(k ),  •••  ,  yq(k  +  L  - 

1)]T  and  nq  =  [nq(k),  ■  ■  ■  ,  nq(k  +  L  —  1)]T,  then  using  the 
matrix- vector  notation  we  can  rewrite  (1)  as 

p 

y,  =  EXPg(M)+n«  (2) 

p=i 

where  Xp  is  an  L  x  M  matrix  which  contains  the  waveforms 
transmitted  from  the  p  th  transmit  element.  Defining  X  = 

[Xi,  •••  ,  XP]  andgg  =  [(g(1’g))T,  •••  ,  (g(P,9))T]T  allows 
(2)  to  be  expressed  as 

y,  =  Xg,  +  n,.  (3) 

By  stacking  the  received  waveforms  across  all  the  Q 
receive  elements,  we  create  y  =  [yiT,  •  •  •  ,  yg]T-  Defining 
X  =  Iq  (g)  X,  we  finally  obtain 

y  =  Xg+  n,  (4) 

where  g  =  [gf,  •••  ,  g^]T  and  n  =  [niT,  •••  ,  nQT]T . 
We  assume  the  target  impulse  response  vector  g  is  a  Gaussian 
random  vector  with  zero  mean  and  full  rank  covariance  matrix 
Eg.  In  particular.  Eg  can  be  diagonalized  through  Eg  = 
UAUff,  where  A  is  a  diagonal  matrix  with  each  diagonal 
entry  given  by  a  real  and  nonnegative  eigenvalue  of  Eg,  and 
U  is  the  corresponding  unitary  eigenvector  matrix.  We  also 
assume  components  of  the  noise  vector  n  to  be  independently 
and  identically  distributed  (i.i.d)  and  complex  Gaussian,  with 
mean  0  and  variance  rr^ . 

The  waveform  design  in  [2]  aims  at  optimizing  two 
performance  measures.  One  measure  is  the  mean-square  error 
resulting  from  estimating  g  with  a  linear  MMSE  estimator, 
which  can  be  expressed  as 

^  =  tr{(a~2XHX  +  i:B-1)-1}.  (5) 

The  other  optimization  metric  is  the  conditional  MI  between 
y  and  g  given  the  knowledge  of  X,  which  can  be  written  as 

1  =  log  [det(a-2EgA'pAf  +  1PQM)]  .  (6) 

Given  the  above  two  performance  metrics,  the  optimal 
design  problem  can  be  formulated  as  to  find  a  waveform 
matrix  X  that  either  maximizes  X  in  (6)  or  minimizes  £  in 
(5),  under  a  constraint  on  the  total  transmit  power  which  is 


expressed  as  tr{XH X}  <  LQVq.  It  turns  out  these  two 
criteria  yield  essentially  the  same  optimum  solution  for  XH X, 
which  is  a  fundamental  quantity  specifying  the  waveform 
design.  Such  equivalence  between  these  two  criteria  is  also  true 
for  an  asymptotic  formulation  from  [2]  which  requires  only  the 
knowledge  of  power  spectral  density  of  the  target.  Further,  the 
optimum  waveform  matrix  X  is  given  by  [2,  Theorem  1] 

X  =  'S>  [{rjlpQM  -  crnA_1)+]1/2  •  VH,  (7) 

where  is  an  LQ  x  PQM  matrix  with  orthonormal  columns 
(i.e.,  semi-unitary),  and  77  is  a  scalar  constant  which  satisfies 

tr{ (jflpQM  -  ct2A-1)+}  =  LQV 0.  (8) 

For  the  asymptotic  formulation,  the  optimal  solution  for  X 
will  be  identical  to  (7),  and  the  only  difference  is  that  U  is 
replaced  by  a  unitary  discrete  Fourier  transform  (DFT)  matrix 
[2,  Theorem  2], 

III.  Problem  Formulation 

Noteworthily,  the  condition  given  in  (7)  is  only  a  nec¬ 
essary  condition  for  X  to  attain  the  optimality,  but  not  a 
sufficient  one.  This  is  because  there  may  exist  some  cases 
where  a  solution  X  satisfies  (7)  but  does  not  have  the  desired 
Kronecker  structure.  For  these  cases,  it  can  be  possible  that  the 
equivalence  between  the  MI  and  the  MMSE  does  not  hold  any 
more.  Therefore,  it  is  not  only  interesting  but  meaningful  to 
seek  a  waveform  solution  or  a  better  approximate  for  X  which 
simultaneously  optimizes  the  performance  measures  of  interest 
and  satisfies  the  structure  constraint.  To  be  more  specific,  we 
would  like  to  explicitly  enforce  the  Kronecker  structure  in 
the  waveform  design.  Towards  this  goal,  we  start  from  the 
optimality  condition  stated  in  (7),  and  formulate  a  new  design 
problem  as  follows. 

Let  us  define  IT  as  a  collection  of  LQ  x  PQM  matrices 
X  that  have  the  desired  Kronecker  structure.  That  is,  for  each 
X  £  II,  there  exists  a  matrix  X  €  CLxPM  which  satisfies 
X  -  1q  0  X.  We  define  a  collection  of  semi-unitary  matrices 
as 

r  =  {'i'  g  clqxpqm  .  =  ipQMy 

r  is  essentially  a  Stiefel  manifold  [5,  p.  301],  which  consists 
of  all  sets  of  PQM  ortho  normal  vectors  in  C,j(L  We  further 
define  a  PQM  x  PQM  matrix 

O  —  [(t?IpQM  —  ^A-1)"1"]1^2  •  Uff.  (9) 

The  right  hand  side  of  (7)  is  then  simplified  as  Thus, 
the  new  waveform  design  problem  can  be  stated  as:  to  find  a 
waveform  matrix  X  £  II  (or  equivalently  a  matrix  X)  that  is 
minimally  distant  from  {\Pf2  :  T }  with  respect  to  the 

Frobenius  norm  or  in  the  sense  of  least  squares.  This  problem, 
if  cast  mathematically,  can  be  the  following: 

Problem  1  :  min  II X  —  % 

X,<5> 

s.t.  x  g  n, 

$er. 


IV.  An  Iterative  Optimization  Algorithm 

Note  in  Problem  1 ,  the  constraint  sets  are  nonconvex. 
To  solve  this  problem,  we  present  an  iterative  optimization 
algorithm  which  features  an  alternating  projection  approach. 
Similar  iterative  algorithms  have  been  suggested  in  [6]  [7], 
and  employed  in  [8]  [9],  for  example.  In  fact,  the  method  of 
alternating  projection,  firstly  proposed  by  John  von  Neumann 
in  1933,  has  a  long  history  and  has  found  applications  to 
many  problems  of  different  disciplines.  It  is  an  approach  for 
finding  the  best  approximation  to  any  given  point  in  a  Hilbert 
space  from  the  intersections  of  a  finite  number  of  subspaces, 
by  alternatively  computing  the  best  approximations  from  the 
individual  subspaces  which  make  up  the  intersection  [10]. 

When  it  comes  to  Problem  1,  the  alternating  projection 
approach  can  start  with  a  matrix  in  one  set,  e.g.,  X  £  II,  from 
which  it  computes  a  \P  £  T  to  minimize  the  Frobenius  norm, 
i.e.,  the  distance  between  the  sets  of  {\Pf2  :  \P  €  T}  and  II. 
Next,  the  algorithm  begins  with  the  obtained  'I'  to  identify  a 
x  gu  to  minimize  the  distance.  Then  the  algorithm  repeats 
this  alternating  process  indefinitely  until  a  certain  stopping  rule 
is  satisfied.  This  process  is  illustrated  in  Fig  1. 


Fig.  1:  An  illustration  of  the  alternating  projection  approach. 

To  this  end,  we  summarize  these  procedures  in  Algo¬ 
rithm  1,  where  k  denotes  the  maximum  number  of  iterations. 
Alternatively,  the  algorithm  can  also  be  halted  when  additional 
iterations  appear  futile  in  further  reducing  the  Frobenius  norm. 
Note  that  the  order  of  these  two  subproblems  can  be  swapped. 
Also,  at  the  initialization  stage,  the  matrix,  either  X  or  \P, 
does  not  need  to  satisfy  the  constraint.  It  should  be  clear  that 
Algorithm  1  never  increases  the  Frobenius  distance  between 
successive  iterations.  In  fact,  it  can  be  proved  that  Algorithm  1 
converges  to  a  global  optimum  in  a  weak  sense.  For  detailed 
results  with  regard  to  its  convergence  properties,  readers  are 
referred  to  [7]. 

To  apply  Algorithm  1 ,  we  need  to  solve  two  minimization 
problems  in  step  3  and  step  4.  These  two  problems  belong  to 
the  general  class  of  matrix  nearness  problems  (see  [11]  for  a 
survey).  We  will  discuss  these  two  problems  in  detail  shortly, 
but  for  notational  simplicity  we  drop  the  iteration  subscript. 
Notice  that  we  will  be  able  to  obtain  closed-form  solutions  for 
both  minimization  problems.  This  is  one  important  appeal  of 
Algorithm  1  as  well  as  this  design  problem. 


Algorithm  1  Alternating  Projection 
Input:  An  initial  ^matrix  X0  G  C  LQxPQM 
Output:  X  and  \P 

l:  i  <—  0; 

2:  while  i  <  k  do 

3:  find  a  matrix  £  T  such  that 

'S',:  =  argmin  \\X,  -  \PI2|||.; 
4:  find  a  matrix  Xi+\  G  II  such  that 

Xi+i  =  arg  min  \\X  - 

Aren 

5:  i  <—  i  +  1; 

6:  end  while 

7:  X  =  Xi+ 1; 

8:  \P  =  \P;. 


A.  Procrustean  Transformation 

In  step  3  of  Algorithm  1,  we  need  to  solve  $£  T  which 
minimizes  \\X  —  \P Tl\\'2F  provided  X  is  predetermined.  Recall 
that  {\Pf2  :  \P  G  T}  consists  of  all  possible  rotations  \PII 
of  SI.  Then  this  problem  is  equivalent  to  finding  a  rotation 
SPST  of  f2  which  best  approximates  X  in  the  sense  of  least 
squares.  In  factor  analysis,  this  type  problem,  particulary  when 
*  is  a  unitary  matrix,  is  usually  referred  to  as  the  problem 
of  finding  a  Procrustean  transformation  of  X  [12,  7.4.8],  or 
simply  an  orthogonal  Procrustes  problem  [13,  Sec.  12.4.1]. 
We  summarize  the  solution  to  this  problem  in  the  following 
theorem. 

Theorem  1:  Let  X  G  CLQxPQM  and  Q  G  CPQMxPQM 

be  given.  Let  flXH  =  VGW,f  be  a  singular  value  decom¬ 
position  (SVD)  of  nXH,  where  V  £  CPQMxPQM  and  W  £ 
C LQxPQM  are  unitary  matrices,  and  G  £  R PQMxPQM 
is  a  diagonal  matrix  with  nonnegative  real  numbers  on  the 
diagonal.  Then  the  solution  to  the  problem 

min  Wx-vnwl  (io) 

is  \P  —  WVH. 

Proof:  We  can  expand  the  least-squares  expression  in 
(10)  as  follows 

\\x-*nfF 

=  tr{(x  -  *n){x  -  vr, i)H} 

=  \\x\\2F-2M{tT{*nxH}}-\M2F 

Since  tr{\prJA’if}  =  tr{JlA’H\P},  to  mimize  || X  —  \PS7||p, 
we  need  to  find  a  matrix  \P  which  solves  the  following 
problem 

max  tr{  JT  XH  (P  }  } . 

Using  flXH  =  VGWff,  we  have 

9?{tr{HA’if'P}}  =  3?{tr{G(WH'PV)}}, 

which  is  maximized  when  Wff\PV  =  I pqm  [12].  This  leads 
to  the  solution  \P  =  WV,f,  and  the  value  of  the  maximum 


is  tr{G}.  This  completes  the  proof. 


given  by 


□ 


B.  Kronecker  Product  Approximation 

In  step  4  of  Algorithm  1,  for  a  given  matrix  we  need  to 
find  a  waveform  matrix  X  which  has  the  required  Kronecker 
structure,  i.e.,  X  =  IQ®X.  For  notational  convenience,  we 
use  X  explicitly  instead  of  X  in  this  part,  and  further  define 
H  =  iFfi.  Thus,  the  problem  can  be  rewritten  as 

min  ||H  —  Iq  <g>  X||!  (11) 

x 

To  solve  the  above  problem,  we  can  reformulate  the  Frobenius 
norm  in  (1 1)  by  appropriately  breaking  matrices  X  and  H  into 
blocks,  a  method  suggested  in  [14].  Thus,  we  obtain 

L  PM 

II H  -  Iq  0  X|j|  =  £  53  ||H(fc,  l )  -  X(fc,  l)  •  IQ\\2F,  (12) 

fe= l  i=i 

where  we  use  the  non-boldface  letter  X(fc,  l)  to  denote  the 
( k ,  /)-th  element  of  matrix  X.  If  using  MatLab  colon  notation, 
H(fc,  l)  can  be  expressed  as 

H (k,  l)  =  H (k  :  L  :  LQ,  l  :  PM  :  PQM). 

In  view  of  (12),  we  can  minimize  ||H  —  Iq  ®  X|||,  simply 
through  a  sequence  of  linear  least  squares  optimizations.  In 
such  a  process,  each  element  of  X  can  be  determined.  We 
summarize  this  result  in  the  following  theorem. 

Theorem  2:  The  matrix  X  £  ^lxpm  c|efplec[  by 

X(k,  l )  =  -Ur{H(fc,  ()},  1  <  k  <  L,  1  <  l  <  PM  (13) 
V 

minimizes  ||H  —  Iq  (g>  X||^. 

Proof:  In  (12),  we  have  expressed  ||H  —  Iq  <g> X|||,  as  a 
sum  of  a  series  of  Frobenius  norms,  and  only  one  term  in  the 
sum  depends  on  X(k,  l)  for  a  given  k  and  l.  For  each  of  these 
terms,  corresponding  to  a  given  k  £  [1  ,L\  and  q  £  [1  ,PM], 
we  have 

||H(M)-X(M)-Iq||f 

-  tr{(H(fc,  l)  -  X(k,  l)  •  IQ)(H(fc,  l )  -  X(k,  l )  ■  I q)h} 

=  ||H(fc,  l)\\2F  -  23?{X*(fc,  l)  ■  tr{H(fc,  ()}}  +  Q\X(k,  l) |2 

which  is  minimized  when 

X(M)  =  ^tr{H(M)}. 

Since  X(fc,  l)  is  specified,  X  is  also  determined.  □ 

V.  Numerical  Results 

In  this  section,  we  present  some  numerical  results  to 
illustrate  the  performance  of  the  waveforms  obtained  through 
the  proposed  algorithm.  Let  us  firstly  describe  the  performance 
benchmarks  under  the  MMSE  and  MI  measures.  We  denote 
the  (i,  *)-th  element  of  the  diagonal  matrix  A  as  A(i,  i).  Thus, 
given  the  optimum  solution  (7)-(8),  the  resulting  value  of 
MMSE  can  be  computed  based  on  [2,  Theorem  1],  and  is 


PQM 


£min  —  ^  ' 


A  (M) 


(A(i,  i)<Jn  21J  1)~*~  +  1  " 

Similarly,  the  resulting  value  of  MI  is  [2] 

PQM 

Xmax=  (log(cr“2A(i,z)r;))+. 


(14) 


(15) 


£min  in  (14)  and  Xmax  in  (15)  represent  the  theoretical  optimal 
values  of  the  objective  functions  (5)  and  (6),  respectively. 
However,  notice  that  when  computing  either  £m;n  or  Imax,  the 
waveform  matrix  X  or  X  is  in  fact  not  used.  Also,  recall  that 
in  the  waveform  design,  the  target’s  second-order  statistics, 
i.e.,  Sg  =  UAU,f.  are  required  a  priori.  Here  we  focus 
on  solutions  of  the  form  discussed  in  [3],  the  asymptotic  case 
where  Sg  is  additionally  assumed  to  have  a  Toeplitz  structure, 
and  can  be  approximated  by  its  associated  circulant  matrix.  In 
this  case,  we  have 


U(M) 


1 

,  exp 

VPQM 


-j2n(k  -  1  )(l  -  1) 

PQM 


for  I  £  [1,  PQM],  Thus  the  diagonal  elements  of  A  are  the 
only  knowledge  we  need  about  the  target’s  statistics. 

Our  numerical  experiments  focus  on  a  MIMO  radar 
scenario  of  P  =  Q  =  2  and  M  =  5,  and  we  fix  the  noise 
power  o2  at  2.  For  simplicity,  but  without  loss  of  generality, 
for  each  single  experiment  we  generate  the  diagonal  elements 
of  A  independently  from  a  uniform  distribution  over  the 
range  [0, 10].  We  can  then  compute  fl  through  (9)  where 
a  water-filling  algorithm  is  employed.  When  implementing 
Algorithm  1,  we  start  from  a  matrix  X{]  which  for  convenience 
is  initialized  through  Xq  =  Iq  <g>  Xo,  with  entries  of  Xo 
given  by  i.i.d.  complex  Gaussian  of  zero  mean  and  unity 
variance.  We  halt  Algorithm  1  when  the  difference  between 
subsequent  values  of  | X  —  \l/n||i?  is  less  than  10-5.  After  the 
waveform  solution  X  is  obtained,  we  calculate  the  resulting 
values  of  MMSE  and  MI  by  simply  plugging  X  into  (5)  and 
(6),  respectively.  We  term  these  results  as  £mjn  and  Imax. 

Our  experiments  indicate  a  close  match  between  the 
theoretical  optimal  performance  (denoted  by  £min  and  Imax), 
and  the  one  attained  by  using  the  resulting  waveforms  (denoted 
by  £min  and  Imax).  However,  due  to  limited  space,  we  are 
unable  to  present  all  these  results  here.  Instead,  we  consider 
another  verification  measure  under  both  the  MMSE  and  MI 
metrics,  -  the  root  mean-square  error  (RMSE).  In  this  context, 
RMSE  is  defined  as  the  square  root  of  the  mean  of  the  squared 
relative  errors  over  all  the  trials,  i.e.. 


RMSE  = 


(16) 


where  N  denotes  the  total  number  of  independent  trials  and 
is  set  to  5000;  e,  denotes  the  relative  error  for  the  i  th  inde¬ 
pendent  trial,  and  is  defined  to  be  the  normalized  difference 
between  iP-m ax )  and  £min  i7"max). 


(a)  RMSE  for  MMSE  versus  L  and  Vo 

Fig.  2:  RMSE  under  the  MMSE  and  MI 


V0 


(b)  RMSE  for  MI  versus  L  and  Vo 
metrics  and  over  N  =  5000  independent  trials 


In  Fig.  2a  and  Fig.  2b,  we  plot  results  under  both  the 
MMSE  and  MI  metrics,  respectively,  for  different  values  of  L 
and  Vo-  It  is  observed  that  as  either  L  or  Pq  increases,  the 
performance  difference  between  the  theoretical  optimal  one 
and  the  one  attained  by  using  X ,  in  terms  of  RMSE,  is  reduced 
accordingly.  Moreover,  in  view  of  the  scale  of  those  results 
depicted  in  Fig.  2a  and  Fig.  2b,  we  find  that  such  performance 
difference,  or  the  RMSE,  seems  almost  negligible.  In  fact,  the 
value  of  RMSE  is  also  very  small  in  each  single  experiment. 
This  clearly  demonstrates  the  effectiveness  of  our  proposed  al¬ 
gorithm,  as  well  as  the  close-to-optimal  feature  of  the  obtained 
waveform  solutions  under  the  MMSE  and  MI  performance 
measures.  Additionally,  it  is  worth  making  a  quick  note  here, 
that  the  close-to-optimal  performance  attained  by  the  resulting 
waveform  solutions  is  also  related  to  the  denseness  of  solutions 
and  the  continuity  of  the  performance  measures. 

VI.  Conclusion 

In  this  paper,  we  revisited  the  MIMO  radar  waveform 
design  problem  originally  investigated  in  [2]  [3],  which  fo¬ 
cuses  on  optimizing  two  criteria:  maximization  of  MI  and 
minimization  of  MMSE.  We  formulated  a  new  waveform 
design  problem,  with  an  aim  to  seek  waveform  solutions  which 
not  only  optimize  the  performance  criteria  of  interest,  but 
have  a  specific  Kronecker  structure.  We  presented  an  itera¬ 
tive  optimization  algorithm  to  find  the  waveform  solutions. 
Such  algorithm  features  the  alternating  projection  method, 
and  consists  of  operations  like  the  Procrustean  transformation 
and  the  Kronecker  product  approximation  to  provide  closed- 
form  solutions  to  some  arising  matrix  nearness  problems.  The 
effectiveness  of  this  proposed  algorithm  has  been  demon¬ 
strated  through  numerical  examples.  These  numerical  results 
also  indicate  that  the  waveform  solutions  obtained  through 
the  proposed  algorithm  can  attain  virtually  indistinguishable 
performance  when  compared  to  that  predicted  in  [2]  [3],  As 
a  result  of  this,  it  is  reasonable  to  conclude  the  equivalence 


between  the  MI  and  MMSE  criteria  in  this  problem  setup, 

which  in  fact  corroborates  the  previous  study  in  [2]  [3]. 
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